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ABSTRACT 

Magnetic flux emergence into the outer layers of the Sun is a fundamental mechanism for releasing energy 
into the chromosphere and the corona. In this paper, we study the emergence of granular-sized flux concentra¬ 
tions and the structuring of the corresponding physical parameters and atmospheric diagnostics in the upper 
photosphere and in the chromosphere. We make use of a realistic 3D MHD simulation of the outer layers of 
the Sun to study the formation of the Ca II 8542 line. We also derive semi-empirical 3D models from non-LTE 
inversions of our observations. These models contain depth-dependent information of the temperature and line- 
of-sight stratification. Our analysis explains the peculiar Ca II 8542 A profiles observed in the flux-emerging 
region. In addition, we derive detailed temperature and velocity maps describing the ascent of magnetic bubbles 
from the photosphere to the chromosphere. The inversions suggest that, in active regions, granular-sized bub¬ 
bles emerge up to the lower-chromosphere where the existing large-scale field hinders their ascent. We report 
hints of heating when the field reaches the chromosphere. 

Subject headings: Sun: chromosphere - Sun: magnetic fields - Radiative transfer - Line: formation - Polariza¬ 
tion 


1. INTRODUCTION 

Currently, one of the foremost questions in solar physics 
concerns the mechanisms that transport and release non- 
thermal energy into the outermost layers of the Sun (e.g. 
Mariska 1992; Klimchuk 2006). The answer to these ques¬ 
tions is evidently connected with understanding the fine struc¬ 
turing observed in the chromosphere, because the conversion 
of non-thermal energy into thermal energy occurs at very 
small spatial scales (e.g.,Tian et al. 2014; Hansteen et al. 
2014; De Pontieu et al. 2014a). Magnetic fields and mag¬ 
netic reconnection play important roles in heating the chromo¬ 
sphere and corona (e.g. Galsgaard & Nordlund 1996; Gudik- 
sen & Nordlund 2005). Reconnection can occur as a result of 
the braiding of the previously existing magnetic field by pho- 
tospheric motions, but can be strongly enhanced when new 
magnetic flux emerges through the photosphere (Galsgaard 
et al. 2007). The emergence of magnetic flux into the outer 
layers of the Sun also plays an important role in the formation 
of sunspots and active regions, and is an integral part of the 
solar cycle (Schmieder et al. 2014). 

During the past years, a number of observational stud¬ 
ies have aimed to describe how small-scale magnetic fields 
emerge in the outer layers of the Sun, and to establish ob¬ 
servational constraints on the shape, field-strength and speed 
of ascent through the photosphere and lower chromosphere 
(Martinez Gonzalez & Bellot Rubio 2009; Schmieder et al. 
2014). In the photosphere, great attention has been given to 
magnetic-flux emergence in the quiet-Sun, because the ubiq¬ 
uity of the latter would translate into large amounts of mag¬ 
netic flux being deposited into the outer layers of the atmo¬ 
sphere (see recent studies by Lites et al. 2008; Beck & Rezaei 
2009; Danilovic et al. 2010; Gosic et al. 2014, and references 
therein). 

In regions with relatively high magnetic activity, and there- 
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fore stronger polarization signatures, observers have focussed 
on the coupling between the photosphere and the lower chro¬ 
mosphere and have studied the connectivity of the emerging 
flux (Guglielmino et al. 2010; Ortiz et al. 2014). However, not 
much work has been done using chromospheric diagnostics, 
probably because non-LTE/non-equilibrium conditions make 
it very difficult to translate the observed intensities into the 
underlying physical state of the atmosphere (see, e.g, Socas- 
Navarro et al. 2000; Leenaarts et al. 2012). 

Realistic radiative numerical MHD simulations of flux 
emergence have, in most cases, been restricted below the up¬ 
per photosphere (e.g., Cheung et al. 2008; Tortosa-Andreu & 
Moreno-Insertis 2009) because non-LTE radiative losses must 
be included to reproduce chromospheric conditions. Using the 
non-LTE method developed by Skartlien (2000), Martmez- 
Sykora et al. (2008, 2009), studied for the first time the emer¬ 
gence of magnetic flux in a simulation including a photo¬ 
sphere, chromosphere and corona. They investigated the chro¬ 
mospheric and coronal response to the emergence of magnetic 
flux, reporting the formation of cold magnetic bubbles in the 
upper photosphere that expel chromospheric oscillations and 
pushes the transition region and the corona to heights much 
greater than 2 Mm. 

A similar scenario is described by Ortiz et al. (2014) (here¬ 
after Paper I), who compare Ca II IR observations with the 
atmospheric parameters of a numerical simulation. They find 
observational evidences of these cold bubbles, which had 
been previously found in simulations by Martinez-Sykora 
et al. (2008, 2009). The ascent of a magnetic bubble from the 
photosphere into the chromosphere was described using very 
high spatial-resolution data in the Fe I 46301/46302 and the 
Ca II 48542 lines. Furthermore, in Paper I the authors also 
presented photospheric LTE inversions from the Fe I dataset, 
and a qualitative study of the Ca II data, including a weak- 
field approximation analysis to infer the strength of the mag¬ 
netic field in the chromosphere. 
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In this paper, we investigate the emergence of small-scale 
magnetic flux from the photosphere into the chromosphere, 
continuing the study initiated by Ortiz et al. (2014). Our anal¬ 
ysis is prompted by our desire to know the structuring of phys¬ 
ical parameters in the upper photosphere and lower/middle 
chromosphere during the emergence of a cold magnetic bub¬ 
ble. We also wish to understand the peculiar shape of the 
Ca II 48542 intensity profiles observed in the flux emergence 
events, with a bump in the red wing that looks like an emission 
peak. The former are obtained by performing non-LTE inver¬ 
sions of high spatial-resolution observations, whereas the lat¬ 
ter is investigated using a realistic 3D MHD simulation of the 
outer layers of the sun. 

In Sect. 2 we describe the processing of the observations, 
the inversions, and the computation of synthetic profiles from 
the 3D MHD simulations. In Sect. 3 we analyze the inverted 
models, and we discuss the formation of the Ca II 8542 line 
during the flux emergence, using the 3D simulation. Our con¬ 
clusions are summarized in Sect. 4. 

2. DATA 

2.1. Observations with the SST 

We analyze the same dataset that was used in Paper I. Our 
observations of AR11024 were acquired on 2009-07-05 start¬ 
ing at 09:48 UT, at coordinates S27°, W12°. The data were 
taken with the Swedish 1-m Solar Telescope (SST, Scharmer 
et al. 2013) and the CRisp Imaging Spectro-Polarimeter 
(CRISP, Scharmer et al. 2008) in the Ca II48542 A line (here¬ 
after 48542). 

We sequentially acquired narrow band images sampling the 
48542 profile at 17 line positions in steps of 100 mA between 
A4 = ±800 mA from line center. An additional point was 
measured at A4 = 2400 mA. 

The data were reduced using the CRISPRED pipeline 
(de la Cruz Rodriguez et al. 2015). The images have been 
restored with multi-frame-multi-object-blind-deconvolution 
techniques (MOMFBD, van Noort et al. 2005) and residual 
seeing distortions within the line scan have been compensated 
using methods described in Henriques (2012). 

It is not possible to reach the continuum of the 48542 line 
with CRISP because the prefilter is not wide enough. There¬ 
fore, to scale our observations relative to the continuum of the 
quiet-Sun spatial average at disk-center, we use a 3D HD sim¬ 
ulation (de la Cruz Rodriguez et al. 2011). The idea is to com¬ 
pute a simulated spatio-temporal average of quiet-Sun spectra 
at yu = 0.87 (where our observations were acquired). Then we 
normalize the average intensity profile by the continuum in¬ 
tensity at /j = 1. Finally, we scaled our observations in such 
a way that the average quiet-Sun spectrum at A4 = 2.4 A has 
the same intensity as the normalised synthetic average profile. 
The same technique was used in Paper I and by Watanabe 
etal. (2012). 

In Paper I we described the peculiar intensity profile of the 
48542 line during the flux emergence event. Figure 1 illus¬ 
trates a time-series of the profile within the region of interest, 
indicated in the corresponding images with a white square. 
The profile has three distinctive features compared with a spa¬ 
tial quiet-Sun average: 

1. The line center is clearly blue-shifted. 

2. Close to line center, the line profile is strongly asym¬ 
metric, with a very sharp gradient on the red side, and a 
much more extended blue side. 
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Fig. 1.— From bottom to top, time evolution of the T8542 spectra, spatially 
averaged within the small white box (7x7 pixels) indicated in the filtergrams 
on the right. The tick-mark separation is 1 " in the images. 

3. At approximately +250 mA, there is an emission fea¬ 
ture that is not similarly visible on the blue wing. 

The degree of asymmetry and the exact strength of the three 
features enumerated above suffer slight changes over time, but 
the overall shape of the profile is persistently similar. In our 
analysis, we will try to understand the formation of the 48542 
line and the atmospheric parameters that are responsible for 
these features. 


2.2. non-LTE calculations 

We use the non-FTE inversion code Nicole to compute in¬ 
versions in the 48542 line (Socas-Navarro et al. 2014). The 
code iteratively modifies the physical parameters of a guessed 
model atmosphere, to reproduce observed spectra (Socas- 
Navarro et al. 2000). Nicole solves the non-FTE problem as 
described in Socas-Navarro & Trujillo Bueno (1997) assum¬ 
ing statistical equilibrium, which is reasonable for Ca II lines 
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TABLE 1 

Number of nodes used in each cycle of our inversions for temperature, 

LINE-OF-SIGHT (l.O.S) VELOCITY AND LONGITUDINAL MAGNETIC FIELD ( Bl ong ). 



Cycle 1 

Cycle 2 

Cycle 3 

Temp 

4 

5 

5 

Vl.o.s 

2 

4 

4 

Blong 

0 

0 

1 


(Wedemeyer-Bohm & Carlsson 2011). Our atom model con¬ 
sists of five bound levels plus an ionization continuum, sim¬ 
ilar to that used by Leenaarts et al. (2009). The atomic data 
used in this study has been extracted from the VALD database 
(Piskunov et al. 1995). 

Nicole assumes Zeeman-induced polarization to compute 
the full-Stokes vector, neglecting scattering polarization and 
the Hanle effect. This assumption works fine in regions with 
relatively strong magnetic-field, but in the quiet-Sun all these 
effects must be considered to model Stokes Q and U observa¬ 
tions (Manso Sainz & Trujillo Bueno 2010; Carlin et al. 2012; 
Carlin & Asensio Ramos 2014). 

The radiative transfer equation is solved using a piece- 
wise cubic Bezier spline interpolant both to calculate the non- 
LTE populations (assuming unpolarized light) and to compute 
the final polarized formal solution (Auer 2003; de la Cruz 
Rodriguez & Piskunov 2013). 

The spatially-averaged quiet-Sun spectrum of the 48542 
line shows a reversed C-shape bisector (Uitenbroek 2006b; 
Pietarila & Harvey 2013; Chae et al. 2013). Only recently, 
Leenaarts et al. (2014) realized that the asymmetry in solar 
observations is produced by the presence of different Ca II 
isotopes, apparently a well-known phenomenon in chemically 
peculiar stars (Castelli & Hubrig 2004; Cowley & Hubrig 
2005; Cowley et al. 2007; Ryabchikova et al. 2008). 

We have included the effect of isotopic splitting in our in¬ 
versions using a simple approximation. First, the non-LTE 
problem is solved with a regular 40 Ca atom. Once the atom 
population densities are known, we compute a final formal so¬ 
lution of the radiative transfer equation, assuming that the ab¬ 
sorption coefficient is a weighted sum of Voigt profiles, each 
of them centered at the rest-wavelength of the line for each 
isotope. In our case, the weight is proportional to the relative 
abundance of each isotope (Carlsson 1986). 

Our observations are under-sampled by almost a factor x2 
in the spectral dimension, a necessary trade-off to keep the 
cadence sufficiently high while observing three spectral lines. 
This situation is less than ideal for inversions, where criti¬ 
cally sampled observations constrain the solution much bet¬ 
ter. Therefore we have run a set of tests to assess the optimal 
combination of nodes that allows us to reproduce the observed 
profiles with the lowest number of degrees of freedom that is 
possible. The inversions are performed in several cycles to im¬ 
prove the convergence (Ruiz Cobo & del Toro Iniesta 1992): 
during the first cycle, a reduced number of nodes is consid¬ 
ered. Once the solution cannot be improved further, a second 
cycle is started from that point with more nodes. Table 1 sum¬ 
marizes the number of nodes used in each of the cycles. 

In our first tests with the number of nodes, we realized that 
placing the nodes equidistantly along the depth-scale is not 
an optimal choice. After some trial and error, we found that 
we obtain much better fits whenever, at least, there is a node 
placed at log T 500 -4.5 both in temperature and in velocity. 

We have modified Nicole to compute the response func¬ 
tions (see definitions in Section 3.2) using centered numerical 


derivates instead of the default derivatives based on one-side 
perturbations of the model parameters. This makes the code 
roughly 33% to 50% slower, but it greatly improves the con¬ 
vergence of the inversion that now reaches more consistently 
similar values of x 2 from different randomizations of the ini¬ 
tial model parameters. Additionally, inspired by van Noort 
(2012), we have removed some inversion noise 1 by inverting 
the data once, then smoothing horizontally the parameters of 
the model at the locations of the nodes, and restarting the in¬ 
version again. This trick only seems to help in pixels where 
the inversion clearly got stuck in a local minimum. 

2.3. Synthetic observation from the 3D simulation 

To study the formation of the Ca II 48542 line we use a 
3D magnetohydrodynamical simulation performed with the 
Bifrost code (Gudiksen et al. 2011), which has also been used 
in Paper I. Synthetic full-Stokes 48542 observations have 
been calculated with Nicole, as described in Sec. 2.2. 

The simulation includes the upper convection-zone, photo¬ 
sphere, chromosphere and corona, covering a physical domain 
that extends vertically from zo = -2.6 Mm to z 1 = 14 Mm 
from the surface of the photosphere. Horizontally, the sim¬ 
ulation covers approximately 24 x 24 Mm. This domain is 
discretized in a grid of 504 x 504 x 496 points, with a hor¬ 
izontal step of 47.6 km and a vertical step of 20 km in the 
radiative zone (-400 < z < 5000 km). In the convection zone 
and corona the vertical distance between grid points becomes 
larger with increasing scale height. 

The calculations are performed with periodic boundary 
conditions in the horizontal plane, and open boundary con¬ 
ditions at the top and at the bottom. At the beginning of the 
simulation, magnetic flux is injected at the bottom boundary. 
The injected field is a flux sheet with no twist oriented along 
the y-axis from vo = 4 Mm to x\ = 16 Mm with strength 
B = 3300 G. This injection continues for 1 hour 45 min¬ 
utes before being turned off. The field rises steadily through 
the convection zone to the photosphere before stalling there 
after about an hour. Eventually the field breaks through the 
photosphere in certain locations and into the chromosphere, 
forming expanding bubbles of magnetic field as described in 
greater detail in paper I. Later stages of this simulation were 
used to study small flares by Archontis & Hansteen (2014). 

We note that compared to the latest runs performed with 
Bifrost, this simulation does not include the effect of colli¬ 
sions between neutrals and ions that can influence the temper¬ 
ature stratification in the chromosphere (Khomenko & Col¬ 
lados 2012; Martinez-Sykora et al. 2012; Martinez-Sykora 
et al. 2015), or the effect of hydrogen non-equilibrium ion¬ 
ization that can increase the amount of free electrons in the 
chromosphere, affecting the opacity of chromospheric lines 
(Leenaarts et al. 2007). 

3. ANALYSIS 

3.1. 3D models from non-LTE inversions 

In this section we discuss the results of the non-LTE inver¬ 
sions of the observed Ca II 48542 Stokes profiles. We will 
focus on the inferred temperature and line-of-sight velocity 
maps. 

Figure 2 illustrates the inferred temperature as a function of 
time at logrsoo = -4.4, -3.0, -1.9, -1.0,0.0. We note that the 

1 Noise in the inferred 2D maps partly from non-smooth convergence of 
adjacent pixels. 
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Fig. 2.— Temporal evolution of the temperature inferred with a non-LTE inversion. From left to right, the panels show consecutive time steps from our time 
series. From bottom to top, the panels illustrate the inferred temperature at iso-log T 500 surfaces in the model. At = 0 corresponds to 10:07:16 UT. 
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Fig. 3. — Temporal evolution of the line-of-sight velocity inferred with a non-LTE inversion. From left to right, the panels show consecutive time steps from our 
time series. From bottom to top, the panels illustrate the inferred temperature at iso-log T500 surfaces in the model. At = 0 corresponds to 10:07:16 UT. 
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Fig. 4.— Temporal evolution of the longitudinal component of the magnetic field inferred with a non-LTE inversion. From left to right, the panels show 
consecutive time steps from our time series. At = 0 corresponds to 10:07:16 UT. 
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panels at log T 500 = 0.0 may be an extrapolation of the model 
in the upper layers because our observations do not reach the 
continuum. Even so, those deep photospheric layers show typ¬ 
ical granulation patterns with hot granules, colder intergran¬ 
ules and hot bright-points (though geometrically the latter are 
probably located in deeper layers than the former, see Steiner 
et al. 2001). 

At log T 500 = -1.0 the granulation pattern is reversed show¬ 
ing hotter intergranular lanes than granules (Cheung et al. 
2007). Here the footpoints of the bubble are clearly visible, 
but otherwise, there are no signatures of the bubble in the 
temperature map. Only at logrsoo ~ -1.9, in the center of 
the bubble above the granule, does the temperature show a 
significant drop compared to the surroundings, especially be¬ 
tween At = 1 min and At = 4 min. Later on, the temperature 
rises again in the area where the bubble appeared. 

At log T 500 = -3.0 there is a dramatic change in the struc¬ 
turing of the temperature, and features that resemble fibrils 
start to be visible. Therefore we associate this panel with the 
lower chromosphere. The magnetic bubble reaches these lay¬ 
ers with a temporal delay of roughly 2 minutes compared to 
log T 500 = -1.9, similarly to what was reported in Paper I 
from the analysis of images at different line positions. 

Above the lower footpoint (in the image) there are signa¬ 
tures of heating, quite localized until At = 4 min when they 
extend around the bubble. We speculate that this heating can 
be produced by three different processes: 

• By magnetic interaction of the footpoints with the ex¬ 
isting magnetic canopy of fibrils. 

• By wave dissipation, assuming that waves are chan¬ 
neled through the footpoints and persistently release en¬ 
ergy into the chromosphere. 

• By Joule heating at the boundary of the magnetic bub¬ 
ble, where steep gradients in the magnetic field drive 
electrical currents. 

Unfortunately, we would need very accurate estimates of the 
magnetic field at these locations to study whether the bright- 
enings are due to reconnection and/or to estimate the loca¬ 
tion and effect of electrical currents in our maps. Therefore, 
this aspect must be clarified in future studies, hopefully with 
higher sensitivity data. 

At log T 500 = -4.4 our inversions show a more homoge¬ 
neous temperature structuring, with no signature of the emerg¬ 
ing field. Unfortunately, this event is right at the end of the 
time-series, so we cannot confirm whether the bubble man¬ 
ages to reach the middle/upper chromosphere: in our dataset 
it does not. The canopy may hinder the ascent of the bubble 
and we find evidences of magnetic interaction between the 
emerging flux and the existing field especially at the bound¬ 
ary of the bubble. 

It is also clear that the apex of the bubble must be located 
between logrsoo = -3.0 and logrsoo = -4.4 in our observa¬ 
tions. We note that the panels in Figure 2 may show reduced 
temperature contrast due to 3D non-LTE effect (see de la Cruz 
Rodriguez et al. 2012) or to straylight in the observations. 

Similarly, 2D maps of inferred velocities are presented in 
Figure 3. Our results in the deepest layers of the inverted 
model seem to be an extrapolation of the velocity field at 
around logrsoo ~ -1.5. We are not surprised because the 
Lorentzian (photospheric) wings are likely to be much more 


sensitive to changes in temperature than to changes in the ve¬ 
locity field, since the intensity gradient is not very steep. Close 
to line center, the profiles carry more complete information of 
the velocity field. 

Our inversion results at logrsoo = -1.9 are very similar 
to those inferred using Fe I lines in Paper I. Gas flows up in 
the center of the bubble and drains down at the sides. These 
downflows coincide with the location of intergranular lanes 
of photospheric granulation. At logrsoo = -3.0, the imprint 
of these downflows in the surroundings has almost vanished, 
but in some locations they are visible at the end of the series, 
when the bubble has reached high enough in the atmosphere. 
Finally, at log T 500 = -4.4 the images are dominated by the 
canopy of fibrils, which appear upflowing along our line-of- 
sight. The entire region seems to be moving upwards with 
speeds of -6 km s -1 . 

We now discuss the inferred longitudinal component of 
the magnetic field, shown in Figure 4. The maps at At = 
0,1,2 min do not show any feature in pixels harbouring the 
bubble. However, from At = 3 min to At = 5 min two patches 
of opposite polarities appear and slowly separate during the 
ascent of the bubble. From At = 6 min, one of the footpoints 
is still visible at the bottom of the image. The inferred lon¬ 
gitudinal field has amplitudes peaking ±140 G which seems 
to lay just above out sensitivity threshold. This inversion has 
been performed with a constant magnetic field because the 
polarization amplitudes are too low to attempt to infer gradi¬ 
ents. 

Initially we were surprised to see these opposite polarity 
patches. In Paper I, we used the far wings and the inner core of 
the 8542 line to estimate the photospheric and chromospheric 
fields, respectively. The results show half-moon shaped foot¬ 
points in the photosphere (with opposite polarities), whereas 
in the chromosphere we detected signals only in a very small 
patch above the footpoints, but unfortunately not in the inte¬ 
rior of the bubble. We suspect this is due to the too restrictive 
wavelength ranges we used to compute the weak-field approx¬ 
imation both at the core and at the wings, skipping those parts 
of the line profile where the bulk of the signals are contained. 

Now we discuss the quality of the fits to the observed pro¬ 
files which are reasonably good, but not perfect. Figure 5 
shows an example from the upper opposite polarity patch, at 
At = 3 min. The general aspect of the Stokes I profile is repro¬ 
duced by the inversion code, but the sharpest features in the 
obsevations appear smoothed in the inversion, an expected re¬ 
sult from the discretization of the nodes in the inversion (de 
la Cruz Rodriguez et al. 2012): the location and separation 
of the nodes are not ideal to catch sharp gradients in the at¬ 
mosphere. The Stokes V profiles are very noisy, but the code 
seems to find a solution that is compatible with the observa¬ 
tions assuming a constant magnetic field. 

We focus now on the strongly tilted blue wing of the pro¬ 
file (see Figure 1). We believe that asymmetry is produced 
in the observations by the presence of fibrils. We note that 
the upper panels in Figure 3 show velocities of the order of 
-6 km s -1 in the fibrils above the bubble, and velocities are 
even higher at lower optical-depths. Those velocities corre¬ 
spond to a Doppler shift of approximately -180 mA from the 
rest wavelength, which is compatible with the location of the 
observed asymmetry in the profile. The asymmetry also quan¬ 
tifies the difference in velocity between the emerging field and 
the fibrils above. We note that in this particular geometry and 
heliocentric angle, the fibrils above the bubble appear upflow- 
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Fig. 5. — Exemplary observed (red-crosses) and best-fit (black-solid) pro¬ 
files from the center of the magnetic bubble at At = 3 min, at the upper 
opposite polarity patch. Stokes / (top) and Stokes V (bottom). 


plotted using arrows in all panels, except in the middle-right 
panel, where the arrows indicate the direction of the magnetic 
field in the plane of our slices (Y — Z). For context, we have 
indicated the location of these slices in Figure 6 using a dotted 
white line. 

In the simulation, the emergence of magnetic-flux produces 
a cold up-flowing bubble, where the plasma flow is mostly 
vertical in the center, and turns almost horizontal on the sides. 
At the boundary with the environment, there seems to be a 
draining effect, and plasma flows down. We note that den¬ 
sity seems to be higher within the emerging region than in 
the surroundings, because the gas is compressed at the upper 
boundary of the bubble. 

To relate the local thermodynamical variables in the simula¬ 
tion to changes in the intensity profile, we have computed the 
contribution function to the intensity and the response func¬ 
tion to temperature, to assess whether the emission feature 
observed in the Ca II 8542 intensity profile can be related to 
temperature enhancements at some height in the atmosphere. 

The contribution function (Cf) is computed using the 
monochromatic source function (S f), and therefore, it en¬ 
codes information from many thermodynamical variables that 
implicitly affect the source function: 

CAz) = S x (z)-e- T ^- d ^~, ( 1 ) 

dz 


ing, but a similar event with downflowing fibrils may show a 
different profile asymmetry. 

3.2. The formation of the T8542 line within an emerging 
bubble 

The advantage of the MHD simulation is that we have all 
the thermodynamical properties in hand, which makes it pos¬ 
sible to study the formation of the line. We have been puzzled 
by the peculiar shape of the line profiles inside the bubble. To 
find an explanation for such profiles, our first attempt was to 
inspect the properties of the 3D simulation. 

Figure 6 shows synthetic images in the Ca II T8542 line at 
different wavelengths, for one time-step of the simulation. In 
the leftmost panel we illustrate the line-core intensity for each 
pixel, and in the consecutive panels, monochromatic images at 
AT = -100, -75,0,75 mA relative to rest wavelength. In the 
first panel, we have effectively removed intensity fluctuations 
due to Doppler motions, and therefore, the dark circular shape 
of the flux emerging region becomes clearer. In the monochro¬ 
matic panels, the intensity also changes due to Doppler mo¬ 
tions, but it is clear that spectra within the emerging region 
are blue-shifted, and therefore, that it is darker at blue wave¬ 
lengths. This is similar to what is observed, as described in 
Paper I. However, the synthetic profiles are significantly nar¬ 
rower than those from real observations, with a much steeper 
Gaussian core than in the observations. Therefore, the imprint 
of Doppler shifts is quite strong in the simulation, and in the 
red wing there is an intensity enhancement from Doppler mo¬ 
tions that is not obviously visible in the observations. 

In the upper and middle rows in Figure 7, we show verti¬ 
cal slices of the corresponding snapshot from the MHD sim¬ 
ulation: in the top row, temperature and line-of-sight veloc¬ 
ity, in the middle row, density 2 and magnetic field strength, 
indicated in black for vertical fields and in red for horizon¬ 
tal fields. For completeness, the velocity field has been over- 

2 normalized to the average density at each horizontal plane. 


where r^iz) is the monochromatic optical depth at height z- 
The response function (Rf) is computed by directly per¬ 
turbing the physical quantities in the model, at each layer, 
and quantifying directly the impact of such a perturbation 
on the synthetic profile (see e.g., Landi Degl’Innocenti & 
Landi Degl’Innocenti 1977, Ruiz Cobo & del Toro Iniesta 
1992, Socas-Navarro & Uitenbroek 2004, Orozco Suarez & 
Del Toro Iniesta 2007): 


Mz,T) 


Sh(z) 

5T(z)6z 


( 2 ) 


where A is the emerging intensity and T is the temperature at 
height z in the model. Fossum & Carlsson (2005) proposed an 
equivalent method to compute the response function, which 
is numerically more stable, with slightly different form than 
Eq. 2, that has been used in our study. We have applied a per¬ 
turbation of ST(z) = ±0.01 T(z). 

In the bottom row of Figure 7, we display the T8542 line, 
computed through the column at the center of the bubble, in¬ 
dicated with a dashed vertical line in the panels above. The 
blue curve corresponds to the temperature stratification (in 
the left panel) and to the velocity stratification (in the right 
panel). The line profile seems to have an emission peak in 
the red wing, similar to what is observed, but contrary to in¬ 
tuition, the temperature stratification does not exhibit clear 
enhancements that could be associated with such peak. Ad¬ 
ditionally, we present the response function to the tempera¬ 
ture and the contribution function for each wavelength (Uiten¬ 
broek 2006a). These two functions quantify how sensitive the 
spectral line is to different layers of the model, however, the 
contribution function does not account for radiative coupling 
of the source function to non-local properties of the atmo¬ 
sphere. 

In our calculations, both methods seem to provide similar 
results: far in the wings, the line is mostly influenced by the 
photosphere, between 0 and 500 km. Close to line center, the 
line is sensitive to higher layers. The blue curve represents 
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Fig. 6. — Synthetic intensity images in the T8542 line. The leftmost panel represents the line minimum intensity at each pixel. All the other panels illustrate, 
from left to right, monochromatic images at AT = -100, -75,0,75 mA, respectively. The vertical dotted line indicates the location of the vertical slices displayed 
in Figure 7. 



Fig. 7. — Vertical slice of the temperature (top-left), vertical velocity (top-right), density (bottom-left) and magnetic-field topology (bottom-right) from a snapshot 
of the simulation. The magnetic field panel illustrates the vertical component in black and the horizontal component along the x-axis in red. The dotted vertical 
line in the center indicates the location where the response and contribution functions have been computed. The velocity field has been overplotted with arrows 
in all panels except that showing the magnetic-field strength, where the direction of the magnetic field in that plane has been represented instead. Bottom-left: 
Response function of the intensity to the temperature. Bottom-right: Contribution function to the intensity. In both cases, the scale has been clipped to enhance 
the upper layers. For clarity, the blue curve in the lowermost panels illustrate the corresponding temperature (in the leftpanel) and velocity profiles (in the right 
panel) for that column in the 3D model. 
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the temperature stratification with height (left panel) and the 
line-of-sight velocity profile (right-panel). Their correspond¬ 
ing scales are indicated in the upper axis of each of the panels. 
Both plots indicate that the core of the line is formed within 
the rising magnetic region, which is denser than the surround¬ 
ings in the simulation. 

The contribution function indicates that the emission peak 
in the red wing originates at around 450 km above the surface, 
coinciding with the base of the bubble. The velocity profile 
also changes there, with a steeper velocity gradient inside the 
bubble. The response function shows a complicated sensitiv¬ 
ity to changes in temperature within a large range of heights, 
making it very hard to identify a specific regime responsible 
for the emission. 

Therefore, we have attacked the problem with a slightly dif¬ 
ferent approach. To test whether the emission is an effect an 
effect of Doppler velocities, we have run a test where the pro¬ 
file from this column is calculated using the original velocity 
stratification, with the velocities set to zero, and with the ve¬ 
locities decreased by a factor x2. The result is illustrated in 
Figure 8 (top panel). The red curve is computed neglecting 
the velocity field, and therefore it is symmetric. The emission 
feature is then visible on both wings of the line, although there 
are no temperature enhancements along the line-of-sight. The 
plot also demonstrates that the emission continues to be vis¬ 
ible when the amplitude of the velocities is half of that in 
the simulations, but the presence of strong Doppler veloci¬ 
ties seem to enhance the emission in the red wing. The latter 
is probably because the velocity gradients (upflowing mate¬ 
rial) reduce the line opacity in the red wing (see discussion by 
Scharmer 1984). 

To understand the origin of the emission, we have inspected 
the shape of the monochromatic source function, and an¬ 
alyzed its dependence on the local temperature profile. In 
Figure 8 (middle panel) we illustrate the temperature pro¬ 
file of the column and the total source function at AT = 
60 mA, which has been converted to radiation temperature. 
The source function decouples from the temperature profile 
at around 350 km from the surface and peaks at 500 km. 
Then it decreases again monotonically. At that wavelength, 
the r y - 1 layer coincides with the bump in the source func¬ 
tion and therefore the emission feature is formed. However, 
this behaviour is not due to an increase in the temperature 
(unlike de la Cruz Rodriguez et al. 2013) or in the density 
profiles, rather the opposite. In the lower panel in Figure 8, 
we show the departure coefficients for the upper and lower 
levels of the 8542 line as a function of height. At z = 500 km, 
the upper level departs from LTE more vigorously than the 
lower level, explaining the peak of the source function at that 
height. 

Uitenbroek (1989) reports that such a behaviour a can occur 
because the Ca II K and the 8542 lines share the upper level 
while their lower levels are collisionally-coupled. Therefore, 
the source function of the 8542 line is set by the K line and 
displays the same temperature sensitivity as 3933 A. We em¬ 
phasize that the source function in Figure 8 decouples from 
the temperature profile exactly at the temperature drop pro¬ 
duced by the magnetic bubble. Given the behaviour observed 
in the departure coefficients and the shape of the source func¬ 
tion, we have no reason to think that this is not the same phe¬ 
nomenon. 

We now focus on the polarization signals in our synthetic 
observations. Figure 9 illustrates some non-consecutive snap- 



Fig. 8.— Top: Synthetic profiles using the original velocity held (solid- 
black), zero velocities (solid-red) and with the velocity held decreased by 
a factor x2 (dotted-blue). Middle: Temperature prohle (solid-black) inside 
the magnetic bubble and the corresponding total source function (red dashed- 
dotted) at AT = +60 mA. The source function has been converted to radiation 
temperature for comparison with the temperature prohle. Bottom: Departure 
coefficients from LTE populations for the upper level (solid) and the lower 
level (dashed-dotted) of the transition. 

shots from the simulation. The left-column displays the min¬ 
imum intensity of the profile at each pixel. In the middle- 
column we have made a color composite: in orange we repre¬ 
sent max( yfi 2^ + Uj) and in blue max(| V^). In this way, we 

can visualize where the field is predominantly inclined. The 
column on the right illustrates Stokes V in the photospheric 
wing, and provides information about the footpoints of our 
emerging region. In the beginning of the time-series, linear 
polarization is dominant in the picture, suggesting that the 
T8542 line is sampling the top of the bubble where the field 
is mostly inclined. As time goes on and the bubble reaches 
higher layers in the atmosphere, the imprint of vertical fields 
becomes stronger on the sides, where the magnetic field con¬ 
nects with the foot-points. This is also consistent with the re¬ 
sults from LTE inversions (using the Fe I A6 301/6302 lines) 
presented in Paper I. Since we are tracking the maximum sig¬ 
nal of the Stokes profiles, which can change wavelength due 
to Doppler motions or to a change in the formation regime, the 
observables depicted in Figure 9 may be sensitive to slightly 
different heights in the model. 

In our observations we did not detect linear polarization 
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within the magnetic bubble and circular polarisation levels 
were just above the detection limit, probably because the mag¬ 
netic field is highly inclined, and the sensitivity of Stokes Q 
and U to the field strength is much lower than for Stokes V. 
Interestingly, we only detected Stokes V signals when the syn¬ 
thetic Stokes I profile shows emission at around Ad = 60 mA 
(e.g. panel 1 of Figure 10 in Paper I). Therefore, we have in¬ 
spected the polarization levels in the Stokes profiles computed 
from the MHD simulation. We note that the longitudinal com¬ 
ponent of the magnetic field inferred from the observations, 
ranges between 100 - 200 G, and these values are quite simi¬ 
lar to those present in the 3D MHD simulation. 

The synthetic Stokes V profiles show polarization levels 
peaking at 0.2% relative to the continuum intensity for the 
linear and circular components. We have convolved the spec¬ 
tra with a theoretical CRISP transmission profile at T8542 and 
most of the signal is washed away, showing peak values one 
order of magnitude lower than the unconvolved ones. How¬ 
ever, the profiles computed from our simulation are signifi¬ 
cantly narrower than the observed profiles (see discussion by 
Leenaarts et al. 2009, de la Cruz Rodriguez et al. 2012, Ru¬ 
bio da Costa et al. 2014). Therefore, our test with the CRISP 
transmission profile is likely to be significantly more pes¬ 
simistic than reality, and it should only be taken as a lower 
limit. Regions with stronger magnetic-field could produce po¬ 
larization levels of « 0.1%, which is our detection limit. Fig¬ 
ure 10 illustrates this experiment. Incidentally, we note that 
the model shows Stokes Q , U and V signals only at the posi¬ 
tion of the emission feature, i.e., the emission brings out po¬ 
larization signals that would otherwise be hidden. This is also 
what happens in our observations. 

Summarizing, in this section we have used a 3D MHD sim¬ 
ulation to study the formation of the T8542 within a magnetic 
bubble emerging into the solar chromosphere. Our analysis re¬ 
veals that the emission component in the red wing of the line 
is consistent with an increase of the source function peaking 
at z = 500 km. The presence of a strong upflow in the bubble 
makes the profile highly asymmetric, hides the emission fea¬ 
ture in the blue-wing and makes the emission in the red-wing 
stronger. 

4. CONCLUSIONS & DISCUSSION 

In this paper we have studied granular-sized flux emergence 
and the formation of the T8542 line within an emerging bub¬ 
ble. To this end, we have performed a multi-pronged analysis 
using a 3D MHD simulation and non-LTE inversions of real 
observations. This study has been driven by two main goals: 
to explain the observed T8542 profiles and to understand the 
structuring of atmospheric parameters within flux-emerging 
regions. 

We have performed non-LTE inversions of our observa¬ 
tional data, obtaining 3D estimates of the temperature and 
the line-of-sight velocity. Our inferred maps show the tem¬ 
poral evolution of the temperature stratification as the bubble 
emerges from the photosphere and protrudes into the lower 
chromosphere. The magnetic bubble leaves a cold imprint 
in the temperature maps, which is first visible in the mid¬ 
dle/upper photosphere and, with a temporal delay, in the lower 
chromosphere. This result is consistent with the analysis per¬ 
formed in Paper I, using the wings of the 8542 line and LTE 
inversions of the observed photospheric lines. Our inversions 
also reveal two patches where the magnetic field has oppo¬ 
site polarity. These patches separate during the ascent of the 
bubble, and they seem to trace the more vertical field of the 



01 2345601 2345601 23456 

x [Mm] x [Mm] x [Mm] 


Fig. 9. — Synthetic observables from our 3D MHD simulation. Left-column : 
intensity image constructed with the minimum intensity of the line at each 
pixel. Middle-column: color composition showing the maximum total linear 
polarization (orange) and the maximum total circular polatization (blue). To 
compose the images, we have scaled the total circular and linear polariza¬ 
tion to the respective maximum values in the selected subfield. Right-column : 
monochromatic Stokes V at AT = -1.83 A from line center. Time increases 
from top to bottom. 

footprints in the upper photosphere, where the bubble is visi¬ 
ble. The inferred longitudinal field peaks at ±140 G. 

Our analysis of the 3D MHD simulation has provided valu¬ 
able information to understand the peculiarities of the spec¬ 
tral profiles that are observed within the emerging region. The 
emission feature that is present in the red wing of the T8542 
line can be explained by the coupling of the T8542 source 
function with the Ca II K line, that can pump electrons to the 
upper level of the transition. Furthermore, the velocity maps 
from non-LTE inversions allow us to relate the asymmetry ob¬ 
served in the blue-wing of the profile with the canopy of fibrils 
above the bubble. 
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Fig. 10.— From top to bottom, Stokes Q, U and V, convolved with a theoret¬ 
ical CRISP transmission profile (red) and unconvolved (black). The profiles 
correspond to the same column from the models as it was used in Figure 7. 


In the 3D MHD simulation, the magnetic region emerged 
without the presence of an organized large scale magnetic 
field in the chromosphere and could therefore reach higher 
layers. In the observations the bubble seems to get stuck once 
it reaches the canopy of fibrils. We believe this is the reason 
why the bubble is visible at line center in our synthetic Ca II 
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8542 observations, but not in real observations (except per¬ 
haps for weak hints of its presence). 

We identify signatures of chromospheric heating in the 
lower chromosphere, but we have not been able to identify the 
heating mechanism that produces these signatures. Also, we 
have not been able to assess whether the magnetic region con¬ 
tinues rising into the upper chromosphere or if it is destroyed 
by interaction with the existing chromospheric magnetic field 
by reconnection. These aspects must be addressed in future 
studies, perhaps including lines that are sensitive to the upper 
chromosphere and transition region (like the Mg II H Sc K 
lines, see e.g., Leenaarts et al. 2013 and the Si IV lines) and to 
the corona. This will be the subject of Paper III of the present 
series (Ortiz et al. in prep.), where we will analyze simultane¬ 
ous observations of flux emergence events with the SST and 
the IRIS space mission (De Pontieu et al. 2014b). 
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